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SOLUTIONS OF TWO HEAT-TRANSFER PROBLEMS WITH APPLICATION 
TO HYPERSONIC CRUISE AIRCRAFT 
Mark D. Ardema 

Office of Advanced Research and Technology 
Mission Analysis Division 
Moffett Field, Calif. 94035 

SUMMARY 


Solutions are obtained for two initial -boundary value problems of one- 
dimensional heat conduction. The problems concern insulation systems for 
liquid-hydrogen-fueled hypersonic aircraft. The solutions are obtained by 
standard analytical techniques and are used to develop procedures for estimat- 
ing the weight of insulation systems for such aircraft. Numerical results are 
presented and compared with results of a finite difference analysis, and 
agreement is found to be excellent. Finally, the solutions are compared with 
steady-state approximations, and it is concluded that these approximations 
provide convenient weight estimating formulas. 

INTRODUCTION 


Studies of loads and weights of liquid-hydrogen-fueled hypersonic air- 
craft are required as part of comparative mission performance analyses. This 
report presents solutions of two heat-transfer problems encountered in esti- 
mating insulation weight of fuselage< thermal protection systems. The analysis, 
however, is general and may be applied to other cases, such as insulated wings 
and fins . 

Previous studies (refs. 1-6) have indicated that hypersonic cruise air- 
craft are likely to be liquid-hydrogen (LH 2 ) fueled and will probably carry a 
major fraction of the fuel in the fuselage. Since the exterior surfaces of 
such aircraft are at high temperatures and their interior structure is at the 
cryogenic temperature of the fuel, their fuselages will require thermal pro- 
tection systems that will significantly increase the gross weight of the 
vehicles. These systems will be required both to prevent excessive boiloff 
of fuel (called the "wet tank" problem in this report) and to limit structural 
temperatures ("dry tank" problem) . 

In the past, the weights of thermal protection systems have been deter- 
mined either from steady-state heat conduction analyses or, in more detailed 
studies, from numerical solution of the transient heat conduction equation. 

It appears that analytical solutions of the heat equation with application to 
hypersonic aircraft thermal protection systems are not available in the liter- 
ature. This report presents such analytic solutions for the wet tank and dry 
tank problems. 

To illustrate the solutions of the two problems, a representative numeri- 
cal example is presented. The results are compared with results obtained with 
a finite difference heat-transfer computer program. The transient solutions 
obtained in the present analysis are also used to investigate the validity of 



steady-state approximations. These approximations are useful for preliminary 
weight estimates because they may be solved for the insulation thicknesses in 
closed form. Both the magnitudes and the sensitivities of these 
approximations are compared with the transient results. 


PROBLEM FORMULATION 


A typical hypersonic aircraft configuration is shown in figure 1 and 


the principal elements of the fuselage 


^7 


Figure 1.- Hypersonic aircraft configuration. 


of such an aircraft are shown schemati- 
cally in figure 2. The fuselage bend- 
ing loads may be carried either by the 
exterior structure (nonintegral tankage) 
or by the tank structure (integral tank- 
age) . Because fuel is used, not all of 
the tank circumference will be in con- 
tact with the LH 2 for the entire flight. 
Consequently, there are two extreme or 
limiting cases to be considered: 



Figure 2.- Fuselage cross section. 


(1) The wet tank case. The tem- 
perature at the wall is held to that of 
the LH 2 throughout the flight. (This 
case corresponds to the bottom of the 
last tank to be emptied.) 

(2) The dry tank case. The heat 
transferred through the insulation is 
absorbed by the tank structure and the 
tank wall temperature is allowed to 
rise accordingly. (This case corre- 
sponds to the top of the first tank to 
be emptied.) 


The insulation thicknesses for each case are given by the solution to the two 
heat-transfer problems associated with these two cases. Vehicle geometry and 
tank sequencing are then considered in estimating the fraction of vehicle 
surface area over which the individual thicknesses apply in the weight 
calculation. 

The major assumptions of the analysis are as follows: (1) heat transfer 

is by conduction only; (2) circumferential heat transfer is negligible com- 
pared with radial; (3) thickness of insulation is small compared with fuselage 
radius; (4) conductivity of all structural elements is infinitely large com- 
pared with conductivity of insulation; (5) insulation is continuous and homo- 
geneous; (6) thermal constants (e.g., conductivity of insulation, specific 
heat of insulation, specific heat of structure)' are independent of position, 
time, and temperature; (7) exterior is exposed to a square temperature pulse 
(or, equivalently, a step pulse). These assumptions imply that initial- 
boundary value problems of the one-dimensional heat equation are to be solved. 
These solutions are obtained by the standard analytical techniques of separa- 
tion of variables and eigenfunction expansions (cf . ref. 7) . 
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For application to hypersonic cruise aircraft, all of the above assump- 
tions are reasonably well satisfied. Since both radiative and convective heat 
transfer at the inner tank surface are relatively small for the tank tempera- 
tures of interest, their effects may be included in the conductivity; thus 
assumption (1) is appropriate. Assumptions (2) and (3) are valid since for 
such vehicles the fuselage diameter is much larger than the structural and the 
insulation thicknesses. For materials usually considered for hypersonic vehi- 
cles, namely metals (cover panels, tank, load bearing structure) and quartz 

fiber type insulations, assumption (4) is 
reasonable. Thus the exterior structure, and 
in the wet tank case, the tank structure, may 
be neglected in the heat-transfer analysis. 
Assumption (5) will be violated in actual 
insulation systems because of heat leaks 
caused by cover panel supports and edge 
effects of finite-dimensional insulation 
blankets. However, these effects may be 
accounted for by using an overall insulation 
conductivity. Although insulation properties 
will typically be a function of temperature 
(fig. 3 ^shows conductivity data for a repre- 
sentative insulation material), representa- 
tive values may be used in order to comply 
with assumption (6) . Later in this paper a 
numerical example shows that solutions 
obtained using a representative value of con- 
ductivity agree reasonably with a variable 
conductivity calculation. Assumption (7) 
requires that the ascent and descent phases 
of the flight take negligible time compared 
with the cruise time. Although this is not 
strictly true, the actual temperature history 
may be approximated reasonably well by a 
square pulse as shown in figure 4. The 
,! actual n pulse shown on this figure was 
obtained from a detailed trajectory computa- 
tion. The n idealized ,! pulse is determined 
by equating the actual pulse area to the 
it should be remarked that the solutions may 


Figure 3.- 


Variation of conductivity with 
temperature. 



Figure 4.- Typical exterior 
temperature pulse. 


idealized pulse area. Finally, 
only be used for single- layer insulations. 


Since the variables of interest (the insulation thicknesses) cannot be 
obtained explicitly, it is necessary to solve for the thicknesses by iteration. 
The thicknesses (weights) are optimized in the sense that the weight of insula- 
tion plus LH 2 boiloff is minimized. The entire analysis is programmed for a 
digital computer. To obtain vehicle thermal protection system weight, items 
such as cover panels, attachments, and "nonoptimum" weight must be added to 
the insulation weight. 

Although application to hypersonic cruise aircraft fuselages is 
exclusively considered in this paper, other possible applications should be 
mentioned. The dry tank solution is directly applicable to insulated struc- 
tures such as wings and fins on cruise vehicles. Its application to vehicles 
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other than cruise vehicles may be limited by assumption (7) above. For 
example, launch vehicles, which have no cruise time, may have exterior temper- 
ature pulses that cannot be approximated by a square pulse. Application to 
reentry vehicles is discussed in reference 8. 


ANALYSIS 


Wet Tank 

The partial differential equation to be solved is 


subject to boundary conditions 


3u , 3 2 u 


u(0,t) = Tg 
u(L,t) = T h 


and the initial condition 


u(x, 0) 



+ 


( 1 ) 

( 2 ) 

(3) 

(4) 



Figure 5.- Wet tank. 


where 



(5) 


This problem is depicted in figure 5. 
Introducing the transformation 


u = 


Th - T S 

v + = x + Tc 


( 6 ) 


into equations (1) - (4) results in the 
following homogeneous version of the 
problem: 


3v 

9t 


, 3 2 v 
k 

3x 2 


(7) 


v(0jt) = 0 
v(L,t) = 0 

v(x,«=I^x + T 0 -^_iSx-T s 


( 8 ) 

(9) 

( 10 ) 
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Setting 


v (x , t) = T(t)$(x) 

in equation (7) results in the separated equations 


dT 

dt 


+ XkT = 0 


( 11 ) 


( 12 ) 


~ = 0 (13) 

, 2 v J 
dx 

Solving equation (13) subject to boundary conditions (8) and (9) gives for the 
eigenvalues and eigenfunctions 

2 

X n = (x) ; n = 1, 2, . . . (14) 

4> n = sin ^ x / X. ~ ^ ; n = 1, 2, . . . (15) 

respectively. Since the solution of equation (12) is 

-\i kt 

T n =e ; n = 1, 2, . . . (16) 

v(x,t) may now be written as 


v(x,t) 



(17) 


n=l 

Since the eigenfunctions are orthogonal for this problem, equations (10) and 
(17) imply 



The temperature distribution is then obtained from equations (6) , (14) , (17) , 

and (18) 
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u(x,t) 



X + 


+ 


oo 


2(T 0 - T S ) 



nir 


-n 2 ir 2 kt 

r 2 


sin 


nnx 

L 


n=i 


(19) 


Since this function is a "strict" solution of the problem it may easily be 
shown that it is also a unique solution. 

The heat transferred through the tank wall to the IH 2 fuel is 


r 


Q = - / f K Mkll dT 
x * 3 x 

'o 


Q = 


K V T s 


T h ) 2 KL(T s - T 0 ) 
+ 


TT k 



(-D 


n 


n‘ 


n=i 


-n 2 ir 2 kt < 


1 - e 


( 20 ) 


An analytical verification of equation (20) is provided by reference 8 which 
gives the solution for the special case of T 0 = Tpj (applicable for re-entry 
vehicles). Setting T 0 = in equation (20) gives 


Q It o =T h 


■ pCL f ~tj ' T S - t h)<i * r t. 
. L p I ,2 PK f 


(- 1 ) 


n 


n 


C .2 2 t 

L P n=l 


, ,2 PK f 

-(mr) z 


1 - e 


C L 2 P 2 


which is equation (2) of reference 8. The heat absorbed by the fuel is 


( 21 ) 


Q “ h fg L B0 


( 22 ) 


where Lgg is the "boiloff thickness" or volume of LH 2 boiled off during 
flight per unit surface area. Hence, combining equations (20) and (22) gives 


Kt f ( T S - T h ) 2KL(T g - T q ) 


J B0 Lh 


£g 


ir 2 kh 


fg 


n=l 


-n 2 Tt 2 ktf 

tULU.-zr- 


n 


( 23 ) 


The insulation weight for the wet tank case can now be computed. The 
weight per unit area is 


6 



W = pL + Ph l BO 


(24) 


where Lgg as a function of L is given by equation (23) . The necessary 
condition for minimum weight is dW/dL = 0 which leads to 


P + PH 


dL 


BO 


dL 


= 0 


V 

PH 


where 


Kt f (T s 


T H> 2K ( T s - T q ) 



(- 1 ) 



2kt, 


-n 2 7r 2 kt^~ l 


n= 1 


= 0 


(-D 


(25) 


n=i 


n 


yj was used. 


Equation (25) must be solved iteratively 


for L. After this has been done, other quantities such as Lgg, Q, and the 
temperature distribution u(x,t) may easily be computed. Computationally it 
was found that all of the infinite series involved converge very rapidly. 

Thus only the first five terms of each series were retained. It is of inter- 
est to note that for fixed materials and temperatures, equation (25) implies 
that for optimum insulation thickness the ratio tf/L 2 is fixed (i.e., the 
optimum insulation thickness is proportional to the square root of cruise 
time) . 


The commonly used steady-state equation may be readily obtained from 
equation (25) by setting T 0 = Tg (from eqs. (1), (2), (3), and (4) this gives 
a linear steady-state temperature distribution) with the result 


PH 


Ktf(T s - 




T 2 

L ss 


(26) 


This may be solved for Lgg in closed form. From equations (20) and (23) the 
steady-state heat transfer and boiloff thicknesses are: 


Qss = 


L 


B0 SS 


K V T s - V 
L ss 

KtfCTg - T h ) 
L SS h fg 


(27) 

(28) 


Equation (26) may be deduced from the transient solution in two other ways. 
First, setting the insulation heat capacity, C, equal to zero (i.e., letting 
k ■+ 00 ) in equation (25) leads directly to equation (26). Secondly, letting 
tf ■+ “ in equation (23) leads to equation (28) which, by virtue of 
equation (24) , leads to equation (26) . 
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Dry Tank 


As shown in figure 6, the initial -boundary value problem to be solved is 




U(L,t) 


0 < t < tj 


Figure 6.- Dry tank. 


9u 

at 



(29) 


u(0,t) = Tg (30) 

f K <•' - A 1 .'.!”' 1 "') - V 

o 

(31) 

T - T 

u(x,0) = — L ° - x + T 0 (32) 


The second boundary condition, equation (31), is a statement that heat trans- 
ferred through the insulation at x = L is equal to the heat absorbed by the 
tank structure. To write this condition in differential form we differentiate 
equation (31) with respect to time and use equation (29) to obtain 


where 


LK_ 


9u(L>t) + k 3 2 u(L, t) _ Q 


9x 


ax" 


K„ = 


B c b l b p b l 


(33) 


Making the transformation 

u = v + T s (34) 

in equations (29), (30), (32), and (33) yields a problem with homogeneous 
boundary conditions 


9v 

at 


a 2 v 


3x 


(35) 


v(0,t) = 0 


avj^ti + k = Q 

B 9x ax 2 


(36) 

(37) 
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v(x,0) = 


(38) 


T - 
H 


x + T, 


We proceed as in the wet tank solution by substituting equation (11) into 
equation (35) to get equations (12) and (13). Applying boundary conditions 
(36) and (37) to equation (13) gives expressions for the eigenvalues and the 
eigenfunctions 

cos^L /F") + /F sin(h /F^ = 0; n = 1, 2, . . . (39) 


4> n = sin(x/F^; n = 1 , 2, . . . (40) 

In this case the eigenvalues cannot be solved for explicitly. Each eigenvalue 
may be shown to lie in the interval 

TT(n - 1) < L /F < tt (n - (41) 

From equations (11), (16), (34), and (40) the solution may be written as 


u(x,t) 



(42) 


and it remains only to find the Fourier coefficients, a . 

The computation of the Fourier coefficients is not straightforward in 
this problem because the eigenfunctions are not orthogonal. This may be 
shown by application of Green's theorem which, in addition, gives a formula 
for the nonorthogonal values : 



1 , 2 , 


\ 

In view of equations (37) , (38) , and (40) , eqaution (43) reduces to 


(43) 



$ . dx 
1 k 


D7y L ’v L)! ii,j 


i, j = 1, 2, . . . 


(44) 
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The nonorthogonality of the eigenfunctions means that the temperature distri- 
bution cannot be obtained in explicit form. This situation, however, does 
not inhibit numerical solution since the problem can be reduced to one of 
matrix algebra. Next introduce the quantities 


I. 

1 


I. . 

il 



v(x,0) A (x)dx ; 


^00 Sh (x)dx ; 


i = 1, 2, . . . 
i, j — 1, 2, ,. 


(45) 


(46) 


where I. and I. . are an infinite-dimensional vector and an infinite- 
l ij 

dimensional matrix, respectively. Using equations (38), (40), and (44) allows 
equations (45) and (46) to be expressed as 


h = Li: (T s - V Si " ( L A ) 

D 

V Si "( L /*r) > i = C 2, ... (47) 

hi ” \ [ L " Lib si " 2 ( L A")] ; 1-1,2.... J 

| (48) 

X ij = ” IF" sin ( L ' /r \) sin ( L A j" ) > i / j ! i, j = 1, 2 , . . . | 

B 


Now evaluate the v(x,t) portion of equation (42) at t = 0, multiply both 
sides by $ m (x), ancl integrate from x = 0 to x = L to get 



v(x,0)$ m 


(x)dx = 



L 

$n( x )$n00 dx 


I 


m 


oo 



a I 
n nm 


(49) 


where equations (45) and (46) were used. Equation (49) gives the Fourier 
coefficients for equation (42) . 


Taking the design point of view for the dry tank weight calculation, it 
is desired to compute the insulation thickness (weight) which limits the 
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structural temperature u(L,t) to a specified value Tg. Since the maximum 
u(L,t) occurs at t £ , equation (42) implies 


Z — A kt ~ . 

a n e n sin ^ L /a^ j (50) 

n=l 

which is to be iteratively solved for L. The numerical procedure for solving 
the dry tank problem is as follows: (1) Guess an L. (2) Solve for the A £ 
from equation (39) with the aid of equation (41) . (3) Solve for the a £ from 
equation (49) with the aid of equations (47) and (48) . (4) Test to see if 

equation (50) is satisfied; if not, repeat the procedure. Although it was 
found that the infinite series involved in the solution converged relatively 
slowly, convergence was entirely satisfactory when the solutions obtained by 
using the first four and the first five terms of these series were averaged. 
This is not surprising since all of the series involved are alternating. Com- 
putationally, the averaging was accomplished by dividing the fifth rows and 
columns of I £ and I £ j by two. 

Unlike the wet tank solution, there is no readily apparent steady-state 
approximation to the transient solution, since the boundary condition at x = L 
makes the problem inherently unsteady. However, a quasi-steady state relation 
may be obtained if it is assumed that the temperature distribution is a linear 
function of x at any time t, that is. 


U ss (x,t) = fj(t) + f 2 (t)x (51) 

This is equivalent to setting C = 0 (k = °°) in equation (29). Since Ugg(x,t) 
must satisfy the boundary and initial conditions (eqs. (30), (31), (32)) with 

T o = T S> 

. Kt 

TI T x fT t u C B L B P B L SS 

= T S ‘ t ( T S " T H^ e ( 52 ) 

b SS 

Since Ugg(Lgg,t^) = Tg, equation (52) leads to 


Kt £ 

CgLgfjgLss 



(53) 


which gives Lgg in closed form. The heat transferred to the structure is 
given by 


Q gs - CgLgpg (Tg - T H ) (54) 

To conclude this section, mention should be made of the relation of these 
two solutions to existing solutions. The solution given in this paper of the 
relatively straightforward wet tank case may also be obtained from the analysis 
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given in section 3.4 of reference 9. This reference also discusses problems 
with the boundary condition of equation (33) , but does not give the solution 
of the dry tank case. However, section 3.13 may be used to obtain the solu- 
tion for the special case t h = V For this case, transformation of equa- 
tion 3.13(9) of reference 9 gives the temperature distribution 


kA_t 


u(x,t) = T s + 2 (T q - T s ) 


n=l 


_ x n u ✓ __ v 

0 n + (L 2 K B 2 /k 2 )]e sin ( /5T" x) 

+ (L 2 K B 2 /k 2 )] + (LKg/k) } 


(55) 


Since it is desired to limit the structural temperature to Tg, this equation 
gives 


-kX„t* 


r S " T B + 2 CT 0 - T S ) 


, -ka 

/ 1 + (L 2 K B 2 /A n k 2 )e 

{ L [ (X^k/LKg) + (LKg/k)] + l( 


= 0 


(56) 


where A n are the roots of equation (39). Equation (56), which is to be 
solved iteratively for L, has application to the insulated structure of 
nonfueled portions of hypersonic cruise vehicles. 


RESULTS 


Comparison With a Finite Difference Solution 

In this section an example is presented to illustrate the solutions 
obtained in the previous section and to provide a comparison with a finite 
difference method. The data for the example was chosen to be representative 
of a hypersonic cruise aircraft travelling at Mach 6 with a cruise time of 
about 1-1/2 hours. The insulation is taken to be quartz fiber, and the tank 
structure, which carries the vehicle loads, is aluminum alloy, weighs about 
3 lb/ft 2 , and is limited to 200° F. The temperature distributions for this 
example for the wet and dry tank cases are shown in figures 7 and 8, respec- 
tively. The flight profile for this case yields: T 0 = 70° F, Tg WET = 952° F, 

T s = 632° F, and T H = -424° F. Note that because of angle of attack 
effects, Tg (which corresponds to the underside of the vehicle) is greater 
than Tg . ^ vera S e insulation thermal properties are those at the following 
temperatures : 

%ET = t H 

T DRY = \ ( T s + T o + t B + t h) 
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(57) 

(58) 




0 FINITE DIFFERENCE METHOD, 

t/t f = I.O, CONSTANT CONDUCTIVITY 



Figure 7.- Wet-tank temperature distributions. Figure 8.- Dry-tank temperature distributions. 

Here T^ is used for since the quantity of importance is heat transfer 

to the LH 2 fuel, and this heat transfer is proportional to the conductivity at 
temperature Tg. For the wet tank case, the optimum insulation thickness L 
was 4.89 in.; the boiloff thickness Lgg was 3.10 in.; the unit weight of 
insulation Wj was 1.48 lb/ft 2 ; the unit weight of boiloff Wgo was 
1.15 lb/ft 2 ; and the heat transferred to the fuel was 222 Btu/ft 2 . For the 
dry tank case, the insulation thickness required to limit the tank temperature 
to 200° F, L, was 2.34 in.; the unit weight of insulation Wj was 0.88 lb/ft 2 ; 
and the heat transferred to the structure was 249 Btu/ft 2 . These values are 
summarized in table 1. Qualitatively, as can be seen from figures 7 and 8 the 
temperature distributions are much as would be expected. As mentioned earlier, 
total vehicle insulation weight may be estimated from a weighted average of 
the weights of the wet and dry tank cases. 


TABLE 1.- COMPARISON OF RESULTS OF PRESENT ANALYSIS 
WITH FINITE DIFFERENCE SOLUTION 



Wet tank 

Dry tank 

L, 

in. 

L B0’ 

in. 

w z , 

lb/ ft 2 

W B0’ 

lb/ft 2 

Q. 

Btu/ft 2 

L, 

in. 

w x , 

lb/ ft 2 

Q> 

Btu/ft 2 

t b> 

°F 



Present analysis 

4.89 

3.10 

1.84 

1.15 

222 

2.34 

0.88 

249 

200* 

Finite difference, 
constant K 

4.89* 

3.17 

1.84 

1.18 

227 

2.34* 

0.88 

248 

198 

Finite difference, 
variable K 

4.89* 

3.38 

1.84 

1.25 

242 

2.34* 

0.88 

240 

178 


* Input values 


A finite-difference heat- conduction program was also applied to the above 
problem. This program solves one-dimensional heat conduction problems by 
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finite difference methods. The insulation thicknesses computed by the ana- 
lytic solutions of the present paper were entered into the finite difference 
program and the resulting temperature distributions were compared with the 
analytically derived ones. The temperature distributions at t = tf for con- 
stant conductivity are indicated by the diamonds on figures 7 and 8. It may 
be seen that these distributions are in excellent agreement with the transient 
ones. Since there is no theoretical difference between these two solutions, 
this agreement indicates a high degree of numerical accuracy. The finite dif- 
ference method gave 227 Btu/ft 2 for heat transferred to the fuel in the wet 
tank case (2.0 percent difference as compared with the analytic solution), 

198° F for the final wall temperature of the dry tank case (as compared with 
200° F) , and 248 Btu/ft 2 for heat transfer to the structure in the dry tank 
case (0.2 percent difference) as shown in table 1. 

Since the finite difference method is capable of solving variable prop- 
erty problems, it may be used to test the validity of using constant conduc- 
tivity. If K is input as a function of temperature (using data from ref. 10), 
the resulting temperature distributions at t = tf are denoted by the squares 
on figures 7 and 8. As shown in table 1, the finite-difference, variable- 
conductivity solution gave 242 Btu/ft 2 for heat transferred to the fuel in the 
wet tank case (9.0 percent difference as compared with the analytical solu- 
tion), 178° F for the final wall temperature of the dry tank case (as compared 
with 200° F) , and 240 Btu/ft 2 for heat transfer to the structure in the dry 
tank case (3.5 percent difference). Although the wet tank constant- 
conductivity calculation using T^et given by equation (57) underestimates 
the heat transfer by 9 percent, it only underestimates insulation plus boiloff 
weight_by about 4 percent. The constant-conductivity dry tank calculation 
with Tory selected from equation (58) gives slightly conservative results. 


Comparison With Steady-State Approximations 


From equations (24), (26), and (28), the steady-state approximation of 
the insulation weight (including boiloff) for the wet tank case is 


%S 


WET 



P H P Kt f *- T S 





(59) 


and using equation (53) the dry tank insulation weight is 

pKt 

Woo = 

SS DRY CgLgpg Zn[(T s - T H )/(T S - T B )] 


(60) 


Since, as mentioned earlier, these equations may be viewed as ignoring the 
heat capacity of the insulation, they will give conservative weights. Equa- 
tions (59) and (60) are convenient for weight estimation because they give the 
weight in closed form. In order to assess the validity of these equations both 
the magnitude of the steady-state weights and the sensitivities of the weights 
to the independent parameters were compared with the weights and sensitivities 
obtained from the transient solutions. Table 2 shows the results of applying 
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TABLE 2.- COMPARISON OF RESULTS OF PRESENT ANALYSIS 
WITH STEADY-STATE APPROXIMATIONS 



Wet tank 

Dry tank 

L, 

in. 

L B0’ 

in. 

Wj, 

lb/ft 2 

W BO’ 

lb/ft 2 

Q> 

Btu/ft 2 

L, 

in. 

w x , 

lb/ft 2 

Q, 

Btu/ft 2 

t b> 

°F 

Present 

analysis 

4.89 

3.10 

1.84 

1.15 

222 

2.64 

0.99 

205 

200* 

Steady-state 

approximation 

4.46 

4.51 

1.67 

1.67 

323 

3.56 

1.34 

205 

200* 


* Input values 


both the transient analyses of the present paper and the steady-state approxi- 
mations to an example case. This example uses a different tank material from 
that used in the example of the previous section. For the wet tank case, the 
steady-state approximation implies that the weight of insulation will equal the 
weight of boiled-off fuel for minimum total weight; whereas the transient solu- 
tion indicates that the insulation weight will be greater than the boiloff 
weight. Since the transient solution gives 2.99 lb/ft 2 for the weight of 
insulation plus boiloff and the steady-state approximation gives 3.34 lb/ft 2 , 
the steady-state approximation overestimates the weight by 11.7 percent. The 
steady-state approximation overestimates the dry tank insulation weight by 
35.7 percent. 


Sensitivities of the weight to four of the independent parameters are 
shown in figure 9. In order to compare sensitivities (slopes), the steady- 
state weights were multiplied by 0.894 and 0.741 for the wet and dry tank 

cases, respectively, to make them 
pass through the same nominal 
point as the transient weights. 

As indicated in figure 9(a), the 
steady-state sensitivities of the 
weight to insulation conductivity 
agree quite well with the tran- 
sient ones; the wet tank case in 
fact agrees virtually exactly. 

It is obvious that this parameter 
is significant for determining 
the weight. Figure 9(b) shows 
the sensitivities to a parameter 
that appears only in the tran- 
sient solutions, the insulation 
heat capacity. The weak influ- 
ence of this parameter indicates 
that the insulation weights will 
be relatively insensitive to 




(a) Insulation conductivity. 


00 


Insulation heat 
capacity. 


Figure 9.- Sensitivity of insulation weight to 
selected parameters. 
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(c) Surface temperature. (d) Maximum structural 

temperature. 

Figure 9.- Concluded. 


vehicle weight. It was mentioned earlier 
parameter, cruise time, of the transient ; 
identical for the wet tank case. 


changes in initial conditions. 

The sensitivity to surface tem- 
perature, shown in figure 9(c), 
shows perhaps the poorest agree- 
ment between transient and steady- 
state solutions. This parameter 
is important since it is a func- 
tion of the highly important 
mission parameter, cruise Mach 
number. An important sensitivity, 
that of maximum structural temper- 
ature, is shown in figure 9(d). 
Since insulation weight decreases 
with increasing structural tem- 
perature and structural weight 
increases, this parameter is 
important for minimization of 
that the sensitivities of a fifth 
nd steady-state solutions is 


It may be concluded from figure 9 that the following equation will give 
a first-order estimate of the average insulation weight (including boiloff) 
per unit area if the values of the independent parameters are "reasonably" 
close to the values of the nominal case: 


W ss =0.894 




pKt f( T S WET - t h) 


0.741pKtf 


fg 


2C b l b p b 


~( Ts dry " t h)/( t s dry " Tb )] 


(61) 


In this equation half the surface area of the vehicle is assumed to be 
governed by the wet tank case and half by the dry tank case. As mentioned 
earlier, thermal protection system weight is obtained by adding such items as 
cover panels, attachments, and "nonoptimum" weights to the insulation weight. 


CONCLUDING REMARKS 


Solutions of two initial-boundary value problems of one-dimensional 
conduction heat transfer have been obtained. The application of these solu- 
tions to hypersonic cruise vehicle thermal protection systems was discussed. 

A numerical example was considered and it was found that the solutions agreed 
well with finite difference solutions. The finite difference solution was 
also used to investigate the validity of the constant conductivity assumption, 
and it was found that this assumption is a reasonable one. Comparison of 
weight sensitivities with those of steady-state approximations showed that 
such approximations give reasonable estimates of the sensitivities, but that 
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these approximations overestimate the weight. An equation for weight 
estimation based on steady-state approximations was presented. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Jan. 27, 1970 
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